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We study the quantum behaviour of a binary mixture of Bose-Einstein condensates (BEC) in 
a double-well potential starting from a two-mode Bose-Hubbard Hamiltonian. We focus on the 
small tunneling amplitude regime and apply perturbation theory up to second order. Analytical 
expressions for the energy eigenvalues and eigenstates are obtained. Then the quantum evolution 
of the number difference of bosons between the two potential wells is fully investigated for two 
different initial conditions: completely localized states and coherent spin states. In the first case 
both the short and the long time dynamics is studied and a rich behaviour is found, ranging from 
small amplitude oscillations and collapses and revivals to coherent tunneling. In the second case 
the short-time scale evolution of number difference is determined and a more irregular dynamics is 
evidenced. Finally, the formation of Schroedinger cat states is considered and shown to affect the 
momentum distribution. 

PACS numbers: 03.75.Lm, 67.85. Fg, 74.50.-|-r 



I. INTRODUCTION 

The experimental discovery of Bose-Einstein condensation [l[ in dilute systems of trapped alkali-metal atoms, such 
as rubidium (Rb), lithium (Li), sodium (Na) and ytterbium (Yb)^ has spurred a renewed interest into the investigation 
of macroscopic quantum phenomena and interference effects, allowing for a deeper understanding of the conceptual 
foundations of quantum mechanics . This fascinating research area has been growing up thanks to the high degree 
of experimental manipulation and control Q . Interference between condensates released in a potential with a barrier 
was first observed in 1997 Q and that paved the way for further investigations on the problem of Bose condensates in 
a double well potential. Then Josephson oscillations have been observed in one dimensional optical potential arrays 
0. A single bosonic Josephson junction was produced for the first time in 2005 with Rb atoms and its dynamics was 
experimentally investigated both within tunneling as well as self-trapping regime @ [s| ■ More recently, mixtures of 
^^Rb and ^'^ Rb atoms have been produced and experimentally investigated [9| as well, whose intraspecies scattering 
lengths could be tunable via magnetic and optical Feshbach resonances. Furthermore the realization of heteronuclear 
mixtures of *^i?6 and ^-^K atoms with tunable interspecies interactions [l^l paved the way to the exploration of double 
species Mott insulators and, in general, of the quantum phase diagram of two species Bose-Hubbard model ill|. The 
interplay between the interspecies and intraspecies scattering produces deep consequences on the properties of the 
condensates, such as the density profile and the collective excitations However, the wide tunability of such 
interactions makes a BEC mixture a very interesting subject of investigation, both from experimental and theoretical 
side as a mean of studying new macroscopic quantum tunneling phenomena as well as the interplay between quantum 
coherence and nonlinearity. Indeed novel and richer behaviours are expected in such a multicomponent BEC. 

On the theoretical side, a bosonic Jose phson j unction with a single species of BEC has been widely investigated by 
means of a two-mode approximation (l^[l5| (l6| , within the classical as well as the quantum regime. In the classical 
regime, characterized by large particle numbers and weak repulsive interactions, the Gross-Pitaevskii equation provides 
a reliable description. Within the two mode approximation it reduces to two generalized Josephson equations which 
describe the time evolution of the relative phase and the population imbalance between the wells [1^1 and differ from 
their superconducting counterpart (l7j by the presence of a nonlinear term which couples the variables. Because of 
such a term, a bosonic Josephson junction exhibit a variety of novel phenomena which range from 7r-oscillations to 
macroscopic quantum self-trapping (MQST) fld\. While the 7r-oscillations, as well the usual Josephosn ones, deal with 
a symmetric oscillation of the condensate about the two wells, the MQST phenomenon is characterized by a broken 
symmetry phase with a population imbalance between the wells. In the quantum regime, characterized by smaller 
values of the particle number and strong interactions, an increasing of phase fiuctuations is observed together with 
the suppression of number fluctuations. Furthermore the time evolution is characterized by phase collapse and revival 
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The quantum behaviour of bosonic Josephson junctions has been deeply investigated by means of the usual 



quantum phase model [l9|[20j[2l| as well as by starting from a two-mode Bose-Hubbard Hamiltonian [22| |23l| [24 
In this context the phase coherence of the junction has been characterized by studying the momentum distribution 
[2^ [2^. The generation and detection of Schroedinger cat states has been investigated as well; indeed the presence 
of such a kind of states reflects in the strong reduction of the momentum-distribution contrast [2^] ^6] . 

More recently such a theoretical analysis has been successfully extended to a binary mixture of BECs in a double 
well potential [23 [1^ [1^ [131 • The semiclassical regime in which the fluctuations around the mean values are small 
has been deeply investigated and found to be described by two coupled Gross-Pitaevskii equations. By means of 
a two-mode approximation such equations can be cast in the form of four coupled nonlinear ordinary differential 
equations for the population imbalance and the relative phase of each species. The solution results in a richer 
tunneling dynamics. In particular, two different MQST states with broken symmetry have been found [29| . where the 
two species localize in the two different wells giving rise to a phase separation or coexist in the same well respectively. 
Indeed, upon a variation of some parameters or initial conditions, the phase-separated MQST states evolve towards 
a symmetry-restoring phase where the two components swap places between the two wells, so avoiding each other. 
Recently, the coherent dynamics of a two species BEG in a double well has been analyzed as well focussing on the 
case where the two species are two hyperfine states of the same alkali metal |31j . 

In this paper we study the quantum behaviour of a binary mixture of Bose-Einstein condensates (BEG) in a 
double-well potential starting from a two-mode Bose-Hubbard Hamiltonian. We analyze in detail the small tunneling 
amplitude regime where number fluctuations are suppressed and a Mott-insulator behaviour is established. We 
perform a perturbative calculation up to second order in the tunneling amplitude and study the stationary states and 
the dynamics of the two species bosonic Josephson junction. Finally, the dynamical generation of Schroedinger cat 
states is investigated starting from an initial coherent spin state and shown to affect the time-dependent population 
imbalance and momentum distribution 25] [2^ . We focus on the contrast in the momentum distribution between 
the two wells and show how it vanishes for a two-component cat state. That could be interesting in view of the 
experimental realization of macroscopic superpositions of quantum states [26] [32] . 

The paper is organized as follows. In Section 2 we introduce our model Hamiltonian within the two-mode ap- 
proximation and define the various parameters. Then we adopt the angular momentum representation and focus on 
the small tunneling amplitude regime. In Section 3 we apply perturbation theory in the tunneling amplitude to our 
Hamiltonian and find analytical expressions for the energy eigenvalues and eigenstates up to second order. Section 
4 and 5 are devoted to the study of the quantum evolution of the number difference of bosons between the two 
wells in correspondence of two different initial conditions: completely localized states and coherent spin states. In 
the first case both the short and the long time dynamics is studied and a rich behaviour is evidenced, ranging from 
small amplitude oscillations and collapses and revivals to coherent tunneling. In the second case the short-time scale 
evolution of number difference is determined and a more irregular dynamics is evidenced, with suppression of the 
dominant frequency when the number of bosons increase. Then, Schroedinger cat states are shown to generate as a 
result of the time-evolution of an initial coherent state when the tunneling between the two wells is suppressed, and 
their influence on the contrast in the momentum distribution is studied. Finally, in Section 6 some conclusions and 
outlooks of this work are presented. 

II. THE MODEL 

A binary mixture of Bose-Einstein condensates [1^ [1^ loaded in a double- well potential is described by the general 
many-body Hamiltonian: 



H = Ha + Hb + Hab (1) 



where 

-2 



(2) 



Hb = Jdl^ l^-^^+v^^, -f i^+Vb (T^) ^b^+lj J dl'd^'i^t (T^) ^+ {-t') Utt (T^ - 1^') (1^') {^) (3) 
are the Hamiltonians for bosons of species a and b respectively and 

Hab= J J d-td-t'^+ {-t) ^+ (7^') Uab - 1^') V-a (l^') i^b (l^) , (4) 
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is the interaction term between bosons of different species. For dilute mixtures one can replace the interaction 
potentials Uaaj C^6h and Uab with the effective contact interactions: 

Uaa (T^ - - gaaS - , U,, (l^ - ^ ) = g,,S (T^ - 7^') , (l^ - 1^') = SatS (l^ - 1^') , (5) 

where gaa — and gi,b — '^^1^°'''' are the intraspecies coupling constants of the species a and b respectively, 

TOa and mt being the atomic masses and Oaa, cibb the s-wave scattering lengths; furthermore gab = ^"'^ is the 
interspecies coupling constant, where iriab = j^'+m reduced mass and Uab is the associated s-wave scattering 

length. In this way the Hamiltonian ©-(II]) can be rewritten as: 



Hab = 9ab J d'tlp'^lp^ tpalpb- (7) 

Here Vi (V') is the double well trapping potential and, in the following, we assume 14 ("^) = Vb (T^) = V (T^); tp^ (T^) , 
ipi (jr^), i = a,b are the bosonic creation and annihilation operators for the two species, which satisfy the commutation 
rules: 

[v^, (T>) , ')] = K (^) ' y^t (^')] - 0, (8) 

(^) , (1^')] = '^^.'5 - > ^ J = a> ^ (9) 

and the normalization conditions: 

'd^ IV-.C^)!' = 7V.; i = a,b, (10) 
Ni, i = a,b being the number of atoms of species a and b respectively. The total number of atoms of the mixture is 

N^Na+Nb. 

Now a weak link between the two wells produces a small energy splitting between the mean-field ground state and 
the first excited state of the double well potential and that allows us to reduce the dimension of the Hilbert space 
of the initial many-body problem. Indeed for low energy excitations and low temperatures it is possible to consider 
only such two states and neglect the contribution from the higher ones, the so called two-mode approximation 
[l5|[l6l| . In this way, by taking into account for each of the two species a and b the mean- field ground states (j)g j 0g 
and the mean- field excited states 0°, (p^, the wave functions ipi, i — a,b, can be rewritten as: 



'fpa = ag(f>'i + QeCpt 

'ipb = bg(t>l + b. 



b , (11) 



2 

where J dPr^ ^\(e) = i = cL,b, and a+, 6+ and a^, b^ {ug, bg and Og, be) are the creation (annihilation) operators 
for a particle of the species a, b in the ground and the excited state respectively. They satisfy the usual bosonic 
commutation relations [0^,0"''] = [^ii^j'^] — ^ij- Furthermore (f>g ^, i = a,b are assumed real for simplicity and such 

that (^{(pg)^ 4>Cj = (^{'Pe)^ 4'g'^ "0, which simplifies the calculations. Let us change the basis and switch to the atom 
number states in such a way that the expectation value of the population of the left and right well can be defined. 
The new annihilation operators are o.l — ^ [o-g + Oe), clr ~ (flg — a^) and 6^ = -1= [bg + b^), = ^ i^g ^ ^e) 
for the species a and b respectively, so that the wave functions ([TT|) become: 

i^a = ^a,. irg + re) + {rg -€) „, 

By substituting Equations (|T2|) into the Hamiltonian ©-([T]), after some algebra we obtain its second quantized version 
within the two-mode approximation: 

^ ^ \ (4afl - aWf ' {a+aL + a+a^) + (a+a^ + a+anf + ^ {b+bR - b+bL^ 
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(13) 



+ -Na {Na - 2) («« ^ + + (7V„\ + iV2^ - Na) + -A^. {E^ + El) + 

where A'i — NiL + Nifi^ i = a,b, is the number of atoms of species a and b respectively, expressed as a sum of numbers 
of atoms in the left and right weU. The parameters are defined as foUows: 



El 



a, b 



(14) 



i — a.b 



(15) 



2 Uap 



(16) 



(17) 



9ab 



(18) 



Aab = 45, 



(19) 



El = 4>i 



i — a,b 



(20) 



k1 , - 2ft' 



i — a,b 



(21) 



2 



e,e 



N,9 



3:6 



a, 6; 



I ^ J. (22) 



In Eq. (1131) . the terms proportional to Ej, i — a,b, describe tunnehng of particles of species a and b from one to the 
other well while the terms proportional to E^^ i = a^b, deal with the local interaction within the two wells and the 
terms proportional to SE'' correspond to additional two-particle processes. Finally the terms proportional to Aab and 
Kfj couple the two species and then various constant terms follow, which we will drop for simplicity. 

In this paper we focus on the small tunneling amplitude regime where number fluctuations are suppressed and a 
Mott-insulator behaviour is established, so it is convenient to introduce the angular momentum representation for the 
species a and b as follows: 



btb 



J' 



1 [bpL 



btb 



J' 



\ {bpR - btb 



(23) 



where the operators Jf , Jf , i — x,y, z, obey to the usual angular momentum algebra and the following relations hold: 



2 V 2 



(24) 
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In particular, the components J] = ^ {Nm — Ni^), i = a,b give the difference of the number of bosons of the species 
i, Nil, and iV^/j, occupying the two minima of the double well potential, i. e. the population imbalances, which are 
experimentally observable quantities. Thus Hamiltonian (fT3|) can be cast in the following form: 

_1_ A ( L^^^ _L 1^^^ i^^^ L-*^^ ^ 

-tl'i^ab'Jz'Jz "Jx-Jx \'^e,g rig^e '^g,g '^e.ej ; l^'J; 

where the constant terms have been dropped. Let us now simplify the notation by introducing the following parame- 
ters: 

hb^El Cb^ASE", Kb = 2^, (26) 



T~j , _ f^ab 1 ab _ ab _ ab 



and rewrite the Hamiltonian as: 



H = ^Aa {J!f - KaJ^ + Ca [J^f + ^Afc [j^f - Kb^ + Cb {J^f + 

+kab.KJl~DabJxJl (27) 

Within the experimental parameters range it is possible to show that d <^ Ai,Ki, i = a,b, and Dab ^ Aab '7'|'29(, 
then in the following we put Ca = Cb = and Dab — 0, which corresponds to neglecting the spatial overlap integrals 
between the localized modes in the two wells. In this way the binary mixture of BECs within two-mode approximation 
maps to two Ising-type spin model in a transverse magnetic field. 

In the following we will focus on the symmetric case A^ = Af, = A and Ka — Kb = K because it allows us to 
perform analytical calculations while capturing many relevant phenomena characterizing the physics of the system. 
So the model Hamiltonian (|25|) becomes: 

H = Ho+Hi, (28) 
Ho = ^A (J^f + Ia {J',f + KabJUl (29) 

Hi = -K {J^ + Jl) , (30) 

where, in the small tunneling amplitude regime, Hj is considered as a perturbation. The total Hamiltonian commutes 
with ( J")^ and (j^) , which leads to the conservation of total angular momentum with quantum numbers ja = 
and jb = ^ respectively. So the whole Hilbert space has finite dimension, equal to {2ja -I- 1) (8i {2jb -I- 1) = {Na + 1) <8) 
{Nb + 1), thus it depends on the number of bosons of the species a and b respectively. The whole basis {ItJq) , \mb)} 
is given by the eigenvectors of J° (J" \ma) — ma l^^a)) and (J^ \mb) = nib |"^h)) with = — and 

llLb — 2 ' 2 ■ 

As a first step we need to diagonalize the unperturbed Hamiltonian (j29p . which can be done by performing the 
following 9 — ^ rotation on the operators J", J^: 

0^=aiJ^-a2Jz n —n — ^ 
—1 , , fli — 02 — —m-, 

J , ^ aiJ'^ + a2Jl ^ 

while an analogous rotation needs to be carried out on J^, entering the perturbation (pO)) . As a result we get: 

H=\{^- Kb) (O') ' + i (A + A,,) (j^) ' - (32) 

which, by defining Ai = A — Aab, A2 = A + Aab, and 0\ — O]. = Jl = = can be cast in the final 

form: 

H = Ai{d]}j^ + A2{jff -2KJl. (33) 

In the following Section we will find analytical expressions for the eigenvalues and the eigenvectors up to second order 
by performing perturbation theory in the tunneling amplitude. 
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III. STATIONARY STATES 

In the present Section we apply second-order perturbation theory to the Hamiltonian of Eq. p3|) in the smah 
tunnehng amplitude limit, which allows us to derive analytical expressions for the stationary states of the system. 
In order to pursue this task let us rewrite Eq. p3p in dimensionless form by assuming ^ as unit of energy: 

H = 2(diy + 2\(^J^y -2kJ^, (34) 

where \ — ^ and fc = then take 

Ai Ai ' 

Ho = 2(diy + 2x(^J^y (35) 

as unperturbed Hamiltonian and 

Hi = -2kJl, (36) 

as a small perturbation term. Here J^, i = x,y,z, obey the usual angular momentum algebra and the following 
relation holds: 



2 ^ 7V2 /iVa 
^ 2 V 2 



1 , (37) 



where: 



= (38) 

In principle, the rotated basis ma)} — { — ^ {ma ~ nii,)'^ , \ni2 — \ {nia + mf,)^} of the unperturbed Hamil- 

tonian ((35|) is given by the eigenvectors of O], = ^ (Ol |mi) = mi |mi)) and (J^ jma) — ma jma)) with mi = 
_j2iziLL^ ^ and ma = —il^tikl^ _ _^ (ja+jb) ^ .^^jj^Qgg corresponding eigenvalues are Eml,m2 = 2 (mi)^ -I- 2 A (ma)^. 

The presence of the operator Ol, which docs not commute with the perturbation term makes the problem of 
finding eigenvalues and eigenvectors of the full Hamiltonian (j34p within perturbation theory much more involved. In 
order to simplify the treatment and carry out analytical calculations while retaining the relevant phenomenology, we 
concentrate on the particular case of a binary mixture where the two species are equally populated, i. e. Na = Ni,, 
and have the same population imbalance between the two wells, i. e. m^ = mb- This situation allows us to describe 
the quantum dynamics of the system in correspondence of the MQST regime, for which we need a completely localized 
initial state. That fixes mi = while ma = rua — — could be an even or odd integer depending on 

even or odd, and leads to the following zero-order eigenvalues: E^"^ j_ = 2A (ma)^ — ^ (ma + rn^f' . Each eigenvalue 

is two-fold degenerate, with the only exception of the ground state for Na even, e'"q^ = 0, which is nondegenerate. 

The two-dimensional subspace of degeneracy is spanned by the states |0,±ma) (where |0,±ma) — ±ma |0, ±ma)) 
and the corresponding zero-order eigenvectors are: 

^o°Lr) = = ^(|0,ma)± |0,-ma)). (39) 

By switching on the perturbation term (l36l) it is possible to show that the degeneration is lifted starting from 
the levels with smaller ma; in general the double degeneracy of the zero-order eigenvalues E^j^± will be lifted at 

the 2ma-th order of perturbation theory [2^. By applying perturbation theory [33| up to order fc^, we obtain the 
following corrected eigenvalues: 

= 2^ (-2) + - 4(^^)2_, ; ^2 ^ 1, 2, (40) 



i^S.=2A+^fM^i±ll±l±MZH±l)); even, (41) 



7 



^J'L = ^ T + - ^ (j2 (J2 + 1) - l) ; ^2 odd, (42) 

where j2 = Furthermore, for A^2 even, the nondegenerate ground state |0,0) belongs to the symmetry class of 
1 0,771 J). The corresponding eigenvectors, up to order k^, are given in the Appendix. 

In the following Sections we use the analytical expressions of energy eigenvectors derived in the Appendix, see Eqs. 
((62)) - (|69p . in order to study the quantum evolution of (t)^ = (i (j^ (r) + ('''))): that is the number difference 
of bosons of species a and b between the two wells of the potential. 



IV. DYNAMICS: COMPLETELY LOCALIZED INITIAL STATES 



In this Section we investigate the quantum evolution of the number difference of bosons of species a and b between 
the two wells assuming a completely localized state as initial condition. That could be interesting in order to elucidate 
the quantum behavior of the system in correspondence of the classical MQST regime and to put in evidence new 
phenomena including quantum coherence in a multicomponent system. In such a case we will study both the short 
and the long time dynamics: as a result a rich behaviour emerges, ranging from small amplitude oscillations and 
collapses and revivals to coherent tunneling. Although such a physics is well known for the single component Bose 
Josephson junction, in our case the dynamics shows that the two species can coexist in the same potential well as if 
there would be an attractive interaction between them. 

As a first step let us recall the general formula which gives the time evolution of the mean value of = i ( J^" + J^*) 



E 



n— m„ n—m, 



± 



0^0n' (ho 



ho,n' ) e V / 



(43) 



where t — -^t is the dimensionless time, the sums are over all the eigenvectors 



being 7r72 = or i. 



and (pn are the projections of the initial state j?/; (0)) on the basis 



N2 

2 ■ 



IV'(o))= E 



(44) 



So it is clear how the knowledge of eigenvalues and eigenvectors is enough in order to study the quantum evolution 



Jt 



of J^, the Bohr frequencies involved, -E'en ~ Eo.n', and the corresponding weights 4'*n4>n' \ ht^.r. 

Let us now study the dynamics of the system when all the bosons of species a and b are initially contained in 
one of the two wells of the potential, say the right one, and then the imbalances of the two species coincide, so that 
NaR = Na, NaL — 0, NbR — Nb, N^L — 0; furthermore the two species are equally populated, i.e. Na = Ni,. That 
implies mi — and 7712 — ^ = in our center of mass rotated basis. The corresponding initial condition is: 



IV (0)) = 



(45) 



In order to investigate the short timescales evolution we need to keep terms up to second order in the tunneling 
amplitude k when we compute the weights in Eq. (j43|) . We find that: 



N2 
2 



2A2 (A^2 - 1)' 



[cos(w^t) - 1] 



(46) 



where the frequency involved is: 



E 



(2) 



E 



(2) 



0,^* 0,(^-1 



2A {N2 - 1) 



A {N2f - 4Af2 



(47) 



At short timescales small amplitude oscillations with frequency w^j around the initial condition ( A^2_r = ^''2 , A^2L = 0) 
are observed and that coincides with a strongly self-trapped regime. 
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In order to investigate the dynamics at longer timescales we have to take into account also the small splittings 
AE N2 ± and AE , N2 \±oi the two higher pairs of quasidegenerate eigenvalues which provide two further frequen- 

cies (see Ref. [23| for the derivation): 

Wo — A£/ ± 



0,^ ~ A"2-i 2«2-2(Ar2_i)! 



fc"2-2 (jV3_i)(jv,_ 2) • (48) 



0,fi^-l) ~ A«2-3 2"2-i(Ar2-3)! 



The whole result is: 



{[nf\r)) = f cos(wor) 



4A2 (7V2 - 1) 



A^2 

— [cos(a;ir) - cos(cjor)] 



+2cos(a;^T) cos(-^t) — cos(c<jit) — cos(cjqt) , (49) 



and, by putting wq = wi = 0, the short timescale dynamics, Eq. (|46p . is recovered. Summarizing, at longer timescales 
the two species bosons are still localized in the initial potential well but the quantum dynamics exhibits collapses and 
complete revivals. Indeed the coefhcient cos (^t), which multiplies the higher frequency term cos(w^r), gives rise 
to the beat, which is responsible for the observed collapses and revivals at timescales fixed by wi, as shown in Fig. 
[TJ Finally, at very large timescales determined by the frequency wo all the bosons tunnel coherently back and forth 
between the two traps; only the first term cos(a;oT) is responsible of such a coherent tunneling, since all harmonic 
functions containing the frequency wi and lu^ are small in amplitude and proportional to , thus they are unable to 
transfer bosons from one trap to the other. 

The tunneling dynamics within macroscopic quantum self-trapping regime described above is analogous to that of 
the TT-mode fixed point obtained by the Gross-Pitaevski approach [23|, where the two species localize in the same 
well despite the repulsive interaction between them. Let us finally note that, despite the explicit dependence on A of 
the frequencies (|T7)) - (|48p . the different physics related to the three time scales described above is simply due to the 
energy splitting introduced by the renormalized tunneling for small Kab- Thus in the case of a mixture of BECs with 
equal population the dynamics remains similar to that of a single component BEG, apart the coexistence of the two 
species in the same well. 

As for the experimental detection of the long timescales phenomena (collapses/revivals and coherent tunneling), 
since the time for their appearance is abruptly increased with A^2, this implies a rapid decrease of the characteristic 
frequencies rendering more difficult the observation of the intermediate and long time behavior in current BEG 
experiments. Indeed pure condensates consisting of 1150 ± 150 atoms of loaded in a double well have been 

recently realized 00 thus rendering the detection of the intermediate time behavior possible. Mixtures with a 
number of atoms ranging from 9 x 10^ and 5 x 10^ [^"^ Rh and '^ifsT Q) to 4 x 10* and 9 x lO'' {^^Rh and [l^) 
have also been recently realized, but in this case very small characteristic frequencies are implied. However, these 
phenomena may be relevant for molecular systems where the number of vibrational excited quanta is small. 

In the next Section we will further investigate the dynamics of the system by assuming as initial state a simple 
coherent state and then study the formation of a particular superposition of such coherent states, the so called 
Schroedinger cat states. 



V. DYNAMICS: COHERENT SPIN INITIAL STATES AND SCHROEDINGER CAT STATES 

In this Section we choose as initial condition a simple coherent spin state [s^l and study the short-time scale 
evolution of number difference; in this way a more complex dynamics will appear. Finally, we study the generation 
of Schroedinger cat states; in particular, we focus on the contrast in the momentum distribution and show how it 
vanishes for a two-component cat state. 

Let us start by considering as initial condition the following coherent spin state [33 |: 

N2/2 I 



\^m-C \ (N2^ WA^2 rytan'-(|)e-*'"=^|0,m2), (50) 



where the coefficient C is: 



C = sin^^/2 ( ^ ) cos^^/2 ( ^ ) e-'(^^/2)0, (51) 
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FIG. 1: (Color onliiio)Tinic evolution of the relative boson number difference between the two traps for different timescales. 
The value of is A; = 0.5 and the boson number is N2 = 10. 



and 9 and are two angles characterizing the superposition. The time evolution of the mean value of up to first 
order in the tunneling amplitude k is given by: 



where 



(3)'"w) = ((«f w) + 



k (aniiO) 



A V 2 



N2 



N2/2-I 

' (f -n)!2(2n+l)- 



[Ci [c0s(We''') — 1] + C2 sin {cOeT) 

N2 - 2n 



An = tan2"+i (^0 [cos (Kr + 0) - cos (<^)] - ^^^2n+i (^"^ 



(52) 



(53) 



with frequencies F„ — — E^^l^^^^ = A (4n + 2). Furthermore the coefficients Ci and C2 are given by: 



^^^ (f + l)!(f -l)! '^°^^^MU-3 



tan'' - - 



tan3 (I) 



*^'''2y tan (I) 



, (54) 



Co = 



tan 



tan I 



^-l) sin(3</,)-r^ + l)sin(0) 



(55) 
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FIG. 2: (Color online)Time evolution of the relative boson number difference between the two traps for different boson numbers. 
The value of fc is fc = 0.1, A = 1.3 and = 10. (black-dashed line), = 15. (red-dotted line), ~ 25. (blue-straight line), 
while 9 — n/2 and — n/i. 



for even, and 



Ci = 



N2I 



No 



(f 



i)!(f -i)! 



■ cos {(f)) 



(56) 



N-2 



■ sin {2cj)) 



tan 



1 



for odd, respectively. Finally, for A^2 even, the zero-order mean value / {J^^ {t) j is given by 



tan(f) 

(0) 



92 



(0) \ No 

(r)\ = -^cos(^) 



sin(6i) 



JV2 



^2! 



tan^ I - 



tan 



(f + l)!(f -1)! 
[cos(cjeT) — 1] + 2 sin (20) sin (wct) > , 



where the dominant frequency uie is equal to uje — — ^_ 

N2 odd is: 



(57) 



(58) 



(-^ -I- 1). The corresponding expression for 



(0) 



No 



(r)} = cos ((?) + - 



2 / tan I 



1 /sin(6i) 



N2 



[cos(weT) — 1] — 2 sin (0) sin (we^) 



(59) 



where uie = E 



(2) 



E' 



= 2k 



N2 (N. 



^ + 1) + i. As one can see, the dominant frequency is gradually suppressed 



with the number of bosons N2 = ^ {Na + Nb) This is clearly seen in Fig. [5] where the boson number difference 
between the two traps is plotted for different A'2 (even) values as a function of the dimensionless time r. One also 
notices a decrease of the oscillation amplitude at increasing N2 ■ 

The effect of A is instead shown in Fig. [3] where the short-time dynamics of the boson number difference is analyzed 
for two values of the interspecies interaction. When A increases the amplitude of the oscillations decreases. The 
detection of the mixture dynamics is thus more favorable for values of A smaller than unity. 
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FIG. 3: (Color online)Time evolution of the relative boson number difference between the two traps for = 10. The value of 
k is k — 0.1 while A — 1.3 (black-dashed line) and A — 0.3(green-straight line), while 8 = tt/2 and (ji = 7r/4. 



A. Cat states 



Let us consider the coherent spin state ([50|) : the expectation value of the Hamiltonian ([34]) on such state is given 
by: 

{ijj iO)\H\i; (0)) ^ 2AnV2 - 2fc ^(^2/2)2 - ^2 cos (60) 

where n = —(N2/2) cos9 and has the maximum value for (p — 0,6 — tt/2. This result also corresponds to the mean-field 
result for the energy. Now, starting from the coherent spin state ([50]) we are interested in looking for Schroedinger cat 
states. Such states are quantum superposition of macroscopic states and their realization has already been suggested 
for a single species Bose-Josephson junction in [23|. Also in the case of a Bose-Josephson junction with binary 
mixtures one might realize cat states from the time-evolution of an initially coherent state following a sudden rise 
of the barrier between the two wells. Thus we consider at time t = 0^ a zero inter-well coupling k, i.e. the time 
evolution is governed by the Hamiltonian Hq in Eq. ([55)1 . For each basis vector |0,m2) of the coherent state (|50p . 
the time-evolution is given by |0, m2)(t) — e^*2'r™2*/^|o, TO2), where T = hn/X is the so-called revival time such that 
li^iT)) = \ip{0)). Considering now the times T/2p, p integer, the time evolution of the coherent state is governed by 
the factor exp(— iTrm^/p) which satisfies the property exp(— i7r(m2 +p)'^/p) = (— 1)'' exp(— iTrm^/p), depending on the 
parity of p. For the choice of even p, a discrete Fourier transform leads to the cat state: 

p-i 

\ij{T/2p)) = J2 Ufce"'=^=/P|e-'2^'=/P^A), (61) 

fc=0 

i.e. a superposition of p coherent states, where Uk = l/p^^~^^ e~''^'"2/Pe'2'^'=™/p. In particular, the cat state 
affects the momentum distribution. This dependence could be important to probe experimentally their existence. In 
particular, when considering the two-component cat state, i.e. for the choice p = 2, one obtains that the contrast in 
the momentum distribution, i.e. the expectation value of on the unperturbed state, vanishes p^. Furthermore, 
the amplitude of the intervals of time in which the contrast is zero increases with increasing 7V2 as clearly shown in 
Fig. H 

It should be noted that despite the close similarity in the behavior of the contrast between the single component 
BJJ and the double one, the mixture will be a better candidate for the creation and detection of cat states. In 
fact their creation time is nh/X and since for repulsive interaction between the two species and A > Aab we get 
A = — [(1 -I- A/Aa6)/(1 — A/ Aab)] > 1, such time can be made short enough to render their detection more favorable. 
For example by fixing the ratio of Rh—^'' Rb interaction to Rh—^^ Rb interaction to be 2.13, a parameter accessible 
in the JILA setup the detection time is twice smaller than the case of a single component BEC. 



VI. CONCLUSIONS AND PERSPECTIVES 



In this paper we investigated the quantum dynamics of a Bose Josephson junction made of a binary mixture of BECs 
loaded in a double well potential within the two-mode approximation. We focused on the small tunneling amplitude 
limit and adopted the angular momentum representation for the Bose-Hubbard dimer Hamiltonian. Perturbation 
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T 

FIG. 4: (Color online) Contrast in the time-evolution of (J^) for 9 — 11/2 and (fi — 7r/4 and even number of bosons. The red line 
is for = 10., the blue one for = 14. and the black one for N2 = 20.. The interval in which the contrast is zero increases 
with increasing 7V2. 



theory up to second order in the tunneling amplitude enabled us perform analytical calculations in the symmetric 
case where Aq = A;, = A and Ka = i^b ~ K . In this way we obtained the energy eigenvalues and eigenstates, whose 
knowledge is mandatory in order to investigate the quantum evolution of the number difference of bosons between 
the two potential wells. In order to study the quantum dynamics more easily and analitycally, we restricted to the 
case in which the two species are equally populated and imposed the condition of equal population imbalance of the 
species a and h between the two wells. We concentrated on the two following initial conditions: completely localized 
states and coherent spin states, and found a rich and complex behaviour, ranging from small amplitude oscillations 
and collapses and revivals to coherent tunneling. Finally, we considered the generation of Schroedinger cat states and 
pointed out their influence on the momentum distribution through the vanishing of the contrast. We showed that the 
creation time can be rendered short enough in the case of a mixture in order to render their detection more favorable. 
That could be crucial in order to build up an experimental protocol to produce and detect cat states within such 
systems. 

We stress that in this work we have chosen to study the symmetric case. This allowed us to obtain analytical 
results, while giving rise to the relevant phenomenology which characterizes the physics of the junction. The general 
case of different couplings between the two bosonic species and / or different populations needs to resort to numerical 
calculations and will be the subject of a future publication (S^- Another interesting issue which deserves further 
investigation is a careful analysis of the quantum manifestations of the self-trapping transition and in general of the 
MQST phenomenon in this more general context. 

The complex dynamics of the generalized Bose Josephson junctions investigated in the present paper could be 
experimentally testable within the current technology. For instance, the JILA group recently Q succeeded in producing 
a mixture of ^^Rh and atoms, whose interactions are widely tunable via Feshbach resonances. In particular it is 
possible to fix the scattering length of as well as the interspecies one and to tune the scattering length of ^^Rh. 
That allows one to explore the parameter space in a wide range and also to realize the symmetric regime A^ = Af, = A. 
Because of the high degree of experimental control, such a setup could be employed to reproduce the phenomenology 
described in this work. 
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Appendix: Order eigenvectors 

The eigenvectors of the full Hamiltonian up to order fc^, are: 



VS) . f 1 - 10,0) + |0, l-> + |^y^iOl±i)Mi±l)^ |„.2.) , ,62) 
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^(2) 
'■0,1- 



1 - 



72A2 



(i2 + l)-2]) |0,l-) + ^x/[i2(i2 + l)-2]|0,2-) 



, J [J2 (J2 + 1)-2][J2 (j2 + l)-6] I , 



(63) 



^(2) 



72A2 



1 - [19h ih + 1) - 2]) |0, 1+) -y MZ1±1) |o, 0) + ^Vb2(i2 + 1)-2]|0, 2+) 



6A 



+ 



/b2 02 + l)-2][j2(j2 + l)-6] 



96A2 



|0,3+), 



(64) 



:(2) 



32A2 



J2 {J2 + 1) - 4 



1± 

0'2 



48A2 



J2 (J2 + 1) - 4 



,15' 

J2 (J2 + 1) - ^ 



5± 
"'2 



^yj2(j2 + l)-jT^- 



J2 (j2 + 1) + 



J2 (J2 + 1) 



(65) 



= Am,\0M) + B+, 0,im2 + lf)+B-^ o,(m2-l)*)+C+, 0, (m2 + 2) 



0,(m2-2)±); m2 5^0,l,^, 



(66) 



where the coefficients are defined as: 



. _ , fc" 4j2 (j2 + 1) (m2) + j2 ij2 + 1) - 4 (m2)^ + 3 (mz)^ 
' ~ 4A2 



(4(m2)' 



1 



(67) 



^± Vi2 {j2 + 1) - m2 (m2 ± 1) 



2A 



(2m2 ± 1) 



(68) 



^± k^ v/i2(j2 + l)-m2(m2±l)Vj^ + l)-(m2±l)(m2±2) 

m2 16A2 



(m2 ± 1) (2m2 ± 1) 



(69) 
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